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It has been shown that there are not only transverse but also longitudinal couplings between microwave 
fields and a superconducting qubit with broken inversion symmetry of the potential energy. Using multiphoton 
processes induced by longitudinal coupling fields and frequency matching conditions, we design a universal 
algorithm to produce arbitrary superpositions of two-mode photon states of microwave fields in two separated 
transmission line resonators, which are coupled to a superconducting qubit. Based on our algorithm, we an¬ 
alyze the generation of evenly-populated states and NOON states. Compared to other proposals with only 
single-photon process, we provide an efficient way to produce entangled microwave states when the interactions 
between superconducting qubits and microwave fields are in the ultrastrong regime. 

PACS numbers: 42.50.Dv, 42.50.Pq, 74.50.+r 


I. INTRODUCTION 

Superconducting transmission line resonators can be used 
as quantum data buses, quantum memories, and single mi¬ 
crowave photon detectors Hi. They usually work in the 
microwave regime and can also be used as quantum nodes in 
so-called quantum networks (l3l01- It is well known that the 
entanglement is one of the most important resources for quan¬ 
tum information processing i], and microwave photons play 
a critical role in quantum state control for solid state quan¬ 
tum devices. Therefore, engineering arbitrarily entangled mi¬ 
crowave photon states is a very fundamental issue for both 
solid state quantum information processing and quantum op¬ 
tics |@] on superconducting quantum chips. 

Usually, nonclassical photon states of a single-mode cavity 
field are generated through the interaction between the cavity 
field and the two-level atom. The methods of generating non- 
classical photon states can be classified into two ways. One 
is to engineer appropriate Hamiltonians in different evolution 
durations by tuning experimental parameters when the target 
state is being generated inni. The other one is to obtain 
the target state via appropriately designed measurements o. 
The former one is deterministic, while the latter one is prob¬ 
abilistic and usually has a low probability to succeed. If the 
nonclassical state is generated using natural atomic systems, 
the latter method is usually more practical since most of pa¬ 
rameters are not possible or not easy to be tuned. However, in 
artificial atomic systems, the former method is more appropri¬ 
ate because system parameters can be artificially controlled. 
For example, superconducting quantum circuits (SQCs) ifl^ 
[l^ provide us a very convenient way to deterministically en¬ 
gineer nonclassical states of a single-mode microwave field by 
varying the system parameters iT HlOl] . 

The method of deterministically generating entangled pho- 
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ton states using atomic systems can be tracked to that of gener¬ 
ating entangled phonon states of two vibrational modes 12^ . 
in a trapped ion interacting with laser fields, by using different 
sideband transitions. However, the number of steps in such a 
method ll^ exponentially depends on the maximum phonon 
numbers. A few proposals were put forward to overcome the 
exponential dependence of the pho non number by introducing 
auxiliary atomic energy levels using phonon number 

dependent interactions ll^ . or employ ing multiphonon tran¬ 
sitions of high phonon numbers These methods have 

successfully reduced the number of steps into quadratic poly¬ 
nomials of the maximum phonon numbers. 

The generation of entangled microwave photon states of 
two modes using superconducting qubit has been studied ll2^ 
1^ . where a classically driven superconducting qubit with 
time-dependent frequency is coupled to two microwave fields 
in two separated cavities. The interaction Hamiltonian be¬ 
tween the superconducting qubit and the cavity fields of two 
modes is described by the Jaynes-Cummings model. There¬ 
fore, there is only single photon transition in each step. How¬ 
ever, the photon-number-dependent Stark effects 11251 - 1^ in¬ 
duced by the qubit-field coupling make it possible to inde¬ 
pendently implement operations for photon states. Thus, the 
number of steps also quadratically depends on the maximum 
photon number. 

It has been shown that the superconducting qubit and the 
cavity field can have both transverse and longitudinal cou¬ 
plings when the inversion symmetry of the qubit potential en¬ 
ergy is broken 1^ . The longitudinal coupling can in¬ 
duce multiphoton transitions ll^ in different sidebands as 
in trapped ions ll^ and thus arbitrary photon states of 
a single-mode cavity field can be more conveniently engi¬ 
neered S. Motivated by studies ll2B - [^ . we study a method 
to generate entangled microwave photon states in two sepa¬ 
rated cavities coupled by a superconducting qubit using mul¬ 
tiphoton transitions. We first show that the longitudinal cou¬ 
plings can induce two-mode multiphoton processes similar to 
those in trapped ions lO . and then study an efficient way to 
generate superposed two-mode photon states. 
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FIG. 1: (Color online) Schematic diagram for a driven qubit (in the 
middle with the blue color), which is coupled to two single-mode 
microwave fields of two separated cavities (in the left with the purple 
color and the right with the red color, respectively). The first cavity 
field has the frequency loi and the second one has the frequency 012 . 
The coupling strength is g\ (^ 2 ) between the qubit and the first (sec¬ 
ond) cavity field. The qubit is driven by a classical field (in the mid¬ 
dle with the black color) with the frequency uj and Rabi frequency 
fl. 

The paper is organized as below. In Sec. [HI an effective 
Hamiltonian, similar to that of trapped ions with two vibra¬ 
tional modes is derived, and then different sideband tran¬ 
sitions are discussed. In Sec. mil a new algorithm is intro¬ 
duced to generate arbitrary superpositions of two-mode pho¬ 
ton states. In Sec. EYl we discuss how to choose parameters 
to obtain a high fidelity of the target state. In Sec.[Vl we nu¬ 
merically study the effects of both imperfect control pulses 
and the environment on the generated target state. In Sec. I VII 
the advantages and experimental feasibility of our method are 
discussed. Finally, we summarize our results in Sec. IVIII 

II. THEORETICAL MODEL AND SIDEBAND 
EXCITATIONS 

A. Basic Hamiltonian 

As schematically shown in Fig.[T] we study a system where 
a superconducting qubit (SQ), modeled as a two level sys¬ 
tem, is coupled to two single-mode microwave fields in two 
separated cavities and driven by a classical field. The system 
Hamiltonian can be given by 

H = Hq + Hr+Hg + Hd. (1) 

Here, Hq and are the free Hamiltonians of the SQ and 
the cavity fields, respectively. Moreover, Hg is the interac¬ 
tion Hamiltonian between the SQ and cavity fields, and Hd is 
the interaction Hamiltonian between the SQ and the classical 
field. In the qubit basis, the qubit Hamiltonian is given by 

Hq = hWq^, (2) 

with CTj; = Ip) (e| + |e) {g\ and = \e) (e| - \g) (p|. The pa¬ 
rameter ujq is the qubit frequency. The kets \g) and |e) denote 
the ground and excited states of the qubit, respectively. 


The free Hamiltonian of two cavity fields is given by 
2 

Hr = '^hujia\au ( 3 ) 

where ai (aj) is the annihilation (creation) operator of the (th 
cavity field with its frequency ug and wi ^ lOi- The interac¬ 
tion Hamiltonian between the qubit and two cavity fields is 

2 

= X! cos6> - sin6>) (aj' +ai), (4) 

1=1 

where gi is the coupling strength between the (th cavity field 
and the qubit, and 0 is a parameter which depends on the in¬ 
version symmetry of the qubit potential energy. 

Similarly, the interaction Hamiltonian between the qubit 
and classical field is given by 

Hd = {dz cos 9 — dx sin 9) cos{ujt + (f)), (5) 

where H is the coupling strength (or Rabi frequency) between 
the qubit and the driving field. The parameters ui and (p are the 
driving frequency and driving phase, respectively. 

In Eqs. (HJi and Q, when the qubit potential energy pos¬ 
sesses inversion symmetry, i.e., cos0 = 0, there are only 
transverse couplings between the qubit and cavity fields 1^ . 
If the rotating wave approximation is further made and there 
is no driving (H = 0), Eq. ([TJ is reduced to extensively stud¬ 
ied Jaynes-Cummings model i^. When the qubit potential 
energy possesses a broken inversion symmetry 1^ . i.e., 
COS0 ^ 0, there are both transverse and longitudinal cou¬ 
plings between the qubit and microwave fields. The bro¬ 
ken inversion symmetry of the qubit potential energy can be 
achieved when the bias charge for the charge qubit or the bias 
flux for the flux qubit is tuned off the optimal point ll^ 1^ . 
But for the phase qubit, the inversion symmetry of the poten¬ 
tial energy is always broken if^[^ . Here, we will study a 
general method and not specify a particular qubit. 

We now change the qubit basis into the current basis of the 
flux qubit or the charge basis of the charge qubit. This is 
equivalent to diagonalizing the operator cos 9 — dx sin 9. 
In the new basis, the Hamiltonian in Eq. O becomes 

H = Hq +Hr+ Hg+ Hd. (6) 

Here, the Hamiltonians Hq, Hg, and Hd are given by 

Hq = hh)x^- + 9ujz^-, ( 7 ) 

2 

Hg='^ hgiaz (a] + ai'^ , (8) 

Hd = Mlaz cos {ujt -F </>), (9) 

with Ox = \g) (e| -F |e) {g\, and Gz = \e) (e| - |p) {g\. Here¬ 
after, the parameters uix = uiq sin 9 and uiz = uiq cos 9 are 
called transverse and longitudinal frequencies of the qubit, re¬ 
spectively. The kets \g) = Ry (—9) |g) and |e) = Ry (—0) |e) 
are persistent current states of the flux qubit or charge states of 
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the charge qubit. Here, Ry {(p) = ex-p{—iipd-y/2) is the rota¬ 
tion operator along the y-axis, with dy = —i |e) (p| -I- i \g) (e|. 
The parameter ^ 0 results in longitudinal couplings be¬ 
tween the qubit and microwave fields in Eq. O- Below, we 
will show that ^ 0 can induce two-mode multiphoton pro¬ 
cesses in the qubit, and then use these multiphton processes to 
generate arbitrary superpositions of two-mode photon states. 


In the interaction picture with the free Hamiltonian Hq = 
+ (f^zCTzl^), Equation (flSl) becomes 


C30 oo 

fc oo oo oo 

— 2 2-^ 2-^ 2-^ ®2 

N— — 00 mi—0 1712—0 
m—o 712—0 

+ H.C., (15) 


B. Multiphoton processes and sideband excitations 


To see how the multiphoton processes can be induced by the 
longitudinal coupling when ui^ ^ 0, we now apply a unitary 
transformation 


D = exp 


.i=\ 


( 10 ) 


to the Hamiltonian in Eq. (|6]). Then, we obtain an effective 
Hamiltonian 


where (f) is the coupling strength between 

the qubit and cavity field with each different transition pro¬ 
cess, and its algebraic form is 


71712712 
^ NlTliTli 


= exp 


2 

NCj + ijJz+ ^ [mi 
1=1 


ni)uJi 


t + iNcj) 


X exp 



Jn {x) 


(- 1 ) 


ni+n2 


mi+ni m2+n2 

'll '/2 


mi\nilm2ln2l 


(16) 


Hes =DHD^ = hujz 


Mlcjz cos (wf -I- (j)) 


( 11 ) 


+ ^ nu:ia\ai + ^ + H.c.' 

It is clear that D is the displacement operator |@] of two-mode 
cavity fields. The displacement quantity is rjiaz/^ for the /th 
cavity field. Hereafter, we will call the picture after the oper¬ 
ator D as the displacement picture. The ratios yi = 2giluji 
are called the Lamb-Dicke parameters in analogy to trapped 
ions Emu. 

To understand the classical-field-assisted multiphoton tran¬ 
sitions of two cavity fields in the qubit, we apply to Eq. (dUia 
time-dependent unitary transformation 


Ud (t) = exp ix-^ sin {ojt -f (p) 


with X = 20./Cj. Then, another effective Hamiltonian 

d 


=UdH,sUl - ihUd^Ul 


— -^'X)zUz 


dt 

2 


ai 


1^1 


H —— [Jnctz^Bn (t)+il.c.], (13) 


N — — 00 

can be derived, with the time-dependent term 
2 


Bn it) = exp 


m - ai^ + iN (Cdt + 4>) 

. 1=1 


( 12 ) 


(14) 


Here, Jn = Jn (x) is the Bessel function of the first kind. 
Equation (fT3l) shows that multiphoton transitions with differ¬ 
ent modes can be controlled by the classical field as in trapped 
ions EH]. 


Equation (fTSl l describes the classical-field-assisted two-mode 
multiphoton processes as in trapped ions EHl . The magnitude 
of J^p2ini depends on x, and p;. We find 


I 71712712 I 

I ^ Nttlitii I 


ktvl 


jlTliTli 


Jn 


71712712 

Jn 


(17) 


where the properties of (t) have been studied 

in Ref. 11^ . The specific expression of J™'"' is given by 


mj.nil \ 


(18) 


X exp {i [Nuj + uiz + {nil — ni) loi] t -f iNcp} . 


Similarly to Eq. (fTTl) . the magnitude of J- 
as 


can be rewritten 




kivl 


jmiTii 


Jn 


(19) 


It is clear that \J7"'‘/Jn\ is independent of the reduced 
driving strength x. Erom Eqs. dnii-illlii, we know that both 
I I can be changed by adjusting x and ry 

in a similar way. By introducing new variables ki = mi — ni, 
we expand Eq. ( fTSl l in the Eock state basis, and then have 


B'lnt —tx 


OO 


oo 


N,ki,k2 — — 00 7li—^l 712—^2 

-H.C.. 


T\rkik2 ('j.'i _ _(1) ^(2) 


ni+ki,ni n2+k2,n2 


( 20 ) 


with = max{0,—A:/} and Q = min{n;,n; + fc;}. Here, 
^ri?n' ~ 7'i\ denotes the ladder operator of the fih cav¬ 

ity field. The time-dependent transition element H^jvCiG 
given by 


W; 


kik2 

Nnin2 


it) = n 


kik2 

N7l\7l2 


exp(iA^i''=t), 


( 21 ) 
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with ni, 712 replaced by Ci, C 2 respectively. The complex tran¬ 
sition amplitude detuning are respectively 

^^^2 _ j\j-Q _j_ _j_ _|_ ]^2U)2- (23) 

The parameter is given by 

I Ip 


with 

_ / 1, fc < 0 

k>0 ’ 

r(fc) (n + fc)! (-1)^ z’- 


(25) 

(26) 


Here, l!^'^ (z) is the generalized Laguerre polynomials. It 
is clear that the classical-field-assisted multiphoton transi¬ 
tions can be derived from Eq. (l20l) using different frequency¬ 
matching conditions. 


C. Time evolution operators 

We now give detailed discussions on how to engineer two¬ 
mode multiphoton processes by tuning the driving field. Let 
us assume that the driving field is tuned to satisfy the resonant 
condition 

A^ ^ = Noj -\- uj^ -f kiuj\ k2^2 — 0- (27) 

Then Eq. dlOl i can be reduced to an effective Hamiltonian 
when unwanted terms are neglected, that is. 


Recall that = max{0, —fc;} and C,i = mm{ni,ni + ki} as 
defined previously. Here the new parameters used in Eq. (l29l l 
are respectively 


^Nnin2 


(f)=COs( 


Qkik2 
^ ^Nnin2 



Qkik2 

‘^Nnin2 


okik2 

‘^Nnin2 


it) = 


— J't>^ 


n "2 sin 


Qfclfc2 
^ ^Nnin2 


jk-ik2 

^Nn\n2 



(30) 
t), (31) 
(32) 


As shown in Eqs. (l28l l and (|29] |. \ki\ photons in the kh res¬ 
onator can be either created if fc; > 0 or annihilated if fc; < 0 
while the qubit is flipped up. Similarly, \ki\ photons in the 
fih resonator can be either created if fc; < 0 or annihilated if 
ki > 0 while the qubit is flipped down. Thus different side¬ 
band excitations can be constructed, depending on the values 
of ki and ^ 2 - 

Because the Hamiltonian derived in Eq. ( fTSl l is similar to 
that of the trapped ions ll^ . the algorithm using two-mode 
multiphonon processes in trapped ions can be directly applied 
into our model, and different superpositions of two-mode pho¬ 
tons can be generated. As a special case, two-mode Lock 
states of high photon numbers can in principle be more ef¬ 
ficiently generated with just two steps as single-mode Lock 
states of high photon numbers I®. However, we here design 
a new algorithm via different sideband transitions of low pho¬ 
ton numbers by tuning the driving field with properly selecting 
the parameters ui^, ^x, i^i, and rji. The detailed discussions of 
parameter selection will be given in Sec.lIVI 


III. ALGORITHM FOR STATE GENERATION 






/ . ^‘‘N(iC2"+"ni+ki,ni"n2+k2,n2 

ni=^l 712=^2 

- H.C.. (28) 


.( 1 ) 


12) 


The time evolution operator governed by the Hamiltonian in 
Eq. (|28] | is given by 


jrk\k2 

Un 


00 00 

= E E^: 

ni=^l 712=^2 


00 00 


E 

E 

ni={i 

712=^2 

oo 

E 

OO 

E 

ni=fl 

712=^2 

OO 

E 

OO 

E 

ni=fl 

712=^2 


kik2 (1) (2) 

N(!^i(^ 2^ 99^ niTLi^ n2n2 


-ik\k2 ( 1 ) (2) 

^ C1 C2 ^ 7 ^ 1 + 1 , Tl 1 ^71 2 + ^2 , 2 


;ikik2^ ^ 


.(1) 


^(2) 

n2,n2-\-k2 


J2) 


-f 


E 


^95^nini^ 112 71,2 ’ 


ni< —fci or 712 < — ^2 


-f E (29) 

ni<fei or 712 <^2 


Let us first study a universal algorithm for generating arbi¬ 
trary two-mode microwave photon states using sideband tran¬ 
sitions with the following four Hamiltonians 
and Here, for the compact of notations, we have used 
k to represent —k with k > 0. Eor instance, is actu¬ 
ally with N = —1, ki = 1, and ^2 = —1. Eor dif¬ 
ferent N, fci, and fe 2 , the interaction Hamiltonian and 

its time evolution operator have already been given in 

Sec. Ill Cl Below, we will first study how to generate the target 
state by choosing pulse durations, frequencies, and phases of 
the driving fields at each generation step with different side¬ 
band excitations, and then we will apply our algorithm to the 
generation of NOON states and discuss particular properties 
of the algorithm. 


A. Universal algorithm for generating arbitrary two-mode 
microwave photon states 

We note that the state generation in our algorithm is studied 
in the displacement picture with the unitary transformation as 
shown in Eq. (fTOl i. The arbitrary quantum states, we expect to 
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be generated, is written as 

1^/) = Cmn2\ni,n2)\g), (33) 

ni+n2<iVmax 

where |ni, 712 ) means that the first and second cavities contain 
ni and 77,2 photons,respectively, and \g) means that the qubit 
is in the ground state. Besides, A^max and mean the 

maximum photon number and the probability amplitude on 
the state |ni, 77-2) \g), respectively. We assume that the system 
is initially in the state 


IV'o) = | 0 , 0 )| 5 ). ( 34 ) 

We suppose the target state |7/>/) can be generated by al¬ 
ternately switching on and off the two-mode transitions 

and With the designed time evolution opera¬ 
tors, the state generation procedure can be represented by. 


iV'/) = n iV'o) 


l'=f 


[u}{tf)A/ {tf)Uf{0) 


ul{h)m^{h)U,{ 0 ) | 0 , 0 )| 5 ), 


( 35 ) 


where Pi, G {11,10,01,00} denotes the transition type for the 
7/th step, and ti, is the time duration for the 7/th step. The time 
evolution operator Ui, (t^) is given by 


U,y {tv) =exp 

* 

-^ tCiaJoi -|- W202®2^ tv 

X exp 

ix— sin {oJvtv + 01/) 



(36) 


As discussed above, the transitions of different types can be 
achieved by changing the frequency w of the driving field, 
which is denoted by ui^, for the 7/th step. The phase of the 
driving field for the 7/th step is denoted by cj)^. We can express 
Eq. drsl) in another equivalent form of iteration. 


= Ul (0) iU) A iU) m , (37) 

with It/jq) and Itjjf) given in Eqs. OTt and ( l33b . respectively. 
The ket is the state after the 7 /th step. We note that the 
subscript / of 1-0/) in Eq. dTSl) denotes the number of the fi¬ 
nal step. Equation dJTl) means that the initial state is restored 
from the target state by a composition of sideband transitions 
with proper time durations, frequencies and phases of driving 
fields. It is a recursion algorithm. 

Without loss of generality, we use the maximum photon 
number A^nax = 2 as an example to show our algorithm. 
The more general case with arbitrary A0nax is given in Ap¬ 
pendix. The detailed steps for generating the target state 


10/) =Co2 |0, 2) \g) + Cn |1,1) \g) + C20 |2,0) jg) 

+C01 |0,1) Iff) + Cio |1,0) Ig) + Coo |0,0) jg) , (38) 


with A^rnax = 2 using our recursion algorithm are described 
as the following four procedures. 



FIG. 2: (Color online) Universal algorithm for generating arbi¬ 
trary two-mode superposition state with the maximal photon num¬ 
ber A^max = 2. The 771 and 772 -axis respectively denote the pho¬ 
ton number of the first and second mode. The two-mode photon 
state is denoted by |77i,n2). The qubit state is represented by the 
q-axis with q = g or e respectively denoting the ground state Ip) 
or excited state |e). The state component |77i,772)|g) is represented 
by a block at the location ( 771 , 772 , 5 ). If a state component is oc¬ 
cupied, we color the corresponding block with red; otherwise, the 
block is left uncolored. The arrows respectively represent the “10 ”, 
“Ol”, “ll”, and “00” transitions with transition types labeled aside 
them. The solid arrow indicates a desired population transfer from 
the starting state to the end state, while the dashed arrow indicates the 
inevitable irrelevant oscillation when the desired population transfer 
is implemented, (a) Schematic diagram for transferring the popula¬ 
tions on states |0, 2)|p), |1, l)|p), and |2, 0)|p) to the state |1, 0)|e). 
This is achieved by consecutively using “Ol”, “lO”, “Ol”, and “lO” 
transitions, (b) Schematic diagram for transferring populations on 
states |0, l)|p) and |l,0)|e) to the state |l,0)|p). This is achieved 
by consecutively using “11” and “00” transitions, (c) Schematic dia¬ 
gram for transferring the population on the state 11, 0) |p) to the state 
|0, 0)|e). This is achieved by using a “lO” transition, (d) Schematic 
diagram for transferring the population on the state |0, 0)|e) to the 
state |0, 0)|p). This is achieved by using a “00” transition. 


Procedure (i). As schematically shown in Eig. EJa), from 
the final state |0/), we first transfer the populations in the 
state space spanned by {|77i, 772 ) 15 ) |77i + 772 = 2} to the state 
|l,0)|e). This procedure consists of four steps as schemati- 
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cally shown in below 

|0,2) |5)^|0,l)|e)^|l,l)|g)^ 

|l,0)|e)^ |2,0)|3). (39) 

In Eq. ( l39b . the transition type and step number is la¬ 
beled respectively above and below the arrow. The ar¬ 
row points to the direction of the population transfer. In 
Step i/, the population transfer is accomplished by properly 
tuning Wj/, and After this procedure, we obtain 

the state Itpf-A) which only has popupaltions in the space 
{|0,0) \g ), |0, 0) |e), |0,1) \g ), |1, 0) \g ), |1, 0) |e)}. We note 
that two additional oscillations 

|l,l)|g)^ |l,0)|e)and|0,l)|5)o |0,0)|e). (40) 

will also occur inevitably when the population transfer from 
the state |0, 2) \g) to the state |0,1) |e) is implemented. But 
they do not cause population leakage outside the original 
space and no extra steps should be taken for them. We call 
these oscillations as irrelevant oscillations. For this proce¬ 
dure, these irrelevant oscillations are schematically shown by 
dashed arrows in Fig.|2|a). The irrelevant oscillation can also 
occur in the following procedures and is shown by dashed ar¬ 
rows. 

Procedure (ii). As schematically shown in Fig.|2|b), start¬ 
ing from the state \ 1 jjf- 4 ), we need to transfer the populations 
in the state space spanned by {|0, l)|g), |1, 0 )|e)} to the state 
|1, 0)|(?). This procedure consists of following two steps 


Thus, we obtain the total step number / = 8 by setting /—8 = 
0. Therefore, the generation of the target state with A^max = 2 
needs 8 steps. 

Our algorithm takes a quadratic number of steps while an 
exponential one is required in Ref. j^ . Fet us now analyze 
the reason. Our algorithm employs four interaction Hamilto¬ 
nians, 7T™, and given in Eq. d^ . However, 

four interaction Hamiltonians “da;(T+ -|-H.c.”, -|-H.c.”, 

“a~ + cr'*'”, and “axOLyCr'^ -f H.c.”are employed in Ref. 12^ . 
The former three interaction Hamiltonians between our algo¬ 
rithm and those in Ref. ll20l] are qualitatively identical since 
they convert the same number of bosons for either mode 
when the two-level system is excited. However, the last ones 
show fundamental difference between our algorithm and that 
in Ref. 1^ . because ours creates one boson (photon) of one 
mode but annihilate one boson (photon) of the other when the 
two-level system is excited. But in Ref. 1^ . one boson for 
both modes can be simultaneously created when the two-level 
system is excited. This difference is critical for us to design 
an algorithm which can keep track of the populations with a 
constant total boson (photon) number. Therefore, there is no 
population leakage outside the original space. However, the 
algorithm in Ref. has population leakage. Obviously, if 
the last interaction in Ref. is changed to ‘"a\.aya'^+H.c’\ 
a theoretically equivalent algorithm to ours can also be de¬ 
veloped. In this sense, our algorithm can be regarded as the 
improved version of that in Ref. . 

B. Calculation of controllable parameters 


|0,l)|g)^|l,0)|e)^|l,0)|g). (41) 

/-4 /-5 

After this procedure, we obtain the state |^/)/_6), which only 
has popupaltions in the space {|0,0) \g ), |0,0) |e), |1, 0) |g)}. 

Procedure (ill). This procedure is similar to Procedure (i). 
As schematically shown in Fig. |2jc), starting from the state 
|'0/_6), here we need to transfer the population on the state 
|1, 0)|(?) to the state |0, 0)|e). This procedure consists of only 
one step as below 


|l,0)|5)^|0,0 )|e). (42) 

/-6 

After this procedure, we obtain the state |'0j_7), which only 
has popupaltions in the space {|0,0) \g ), |0,0) |e)}. 

Procedure (iv). This procedure is similar to the Procedure 
(ii). As schematically shown in Fig. |2|d), starting from the 
state \tjjf-j), we need to transfer the population on the state 
|0, 0)|e) to the state |0, 0)|5). This procedure only consists of 
one step as below 

|0,0 )|e)^|0 ,0 )|5). (43) 

We thus obtain the state |'0/-8) = |0, 0)|g). 

Therefore, the target |'0j) can be generated from the initial 
state I'i/’Z-s) using inverse processes from the Procedure (iv) 
to the Procedure (i). We note {ipf-s) = IV'o) = |0, 0)|(7). 


Fet us now study how to choose the pulse duration ti,, the 
frequency ui^ and phase ip^, of the driving field to generate a 
target state in the i^th step for different types of transitions. 

We suppose that the population transfer is taken as follow¬ 
ing 


|ni,n 2 ) \g) \ni + ki,n 2 + k 2 ) \e ), (44) 

V 

in the I'th step, where the transition type Pi, = ki k 2 should 
be switched on based on the previous discussions. Thus the 
driving frequency is taken as 


COi, = UJz + kiLUi + ^2^2, (45) 

from the resonant condition in Eq. dZTl i. By introducing the 
notations 


CiiL = (si IV'.'), (46) 

= («i:«2| (e| ( 47 ) 

then from Eq. dTTl) . we need to solve the equation, 

(ni,n2| (si Ul (0) (f,) tJ, (U) IV'.) = 0. (48) 


We thus have the explicit solution for the pulse duration as 


tu = 


Qfclfc2 

1C1C2 


arctan 


a 




ni,n2 


D, 


ni+fei ,712 + ^2 


(49) 
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The phase of the driving field is determined by 


arg 


a 


(•') 


ni + fci ,n2+k2 ^ 


D 


=x sin {ojyti, + (jiv) 

(50) 


Here, the notation is the value of for the z^th 

step, which is given in Eq. (|32]i and depends on (f)^. Still recall 
0 = min{ni,ni + A:;} with / = 1 , 2 . 

Similarly, if the population transfer is taken as 

|ni + fci,n 2 + ^ 2 ) |e) ^ |ni,n 2 ) \g ). (51) 


in the v\h step. The explicit solution for is then 


tu = 


q/ci/C2 

“ 1 C 1 C 2 


arctan 


L». 


Til +^1 ,^2+^2 


a 




ni ,712 


(52) 


and the phase of the driving field is determined by 


arg 


a 




D 


ni-\-ki,7i2+k2 / 


=x sin (wyfj, + 


^fcife2i/_|—niod27r. (53) 
1C1C2 2 


According to the target state, the time duration, frequency 
and phase of the driving field for each step can be calculated 
using above equations. For example, if the “00” transition is 
used in the 3rd step, then we use Eq. (l52]) and Eq. (l53]) to 
obtain and (^3 by setting v = 3. 


C. Application to NOON states 

As an example, we now apply our algorithm to the genera¬ 
tion of the NOON state, i.e., the target state is 

lijf) = ^ (|iV„,ax,0) \g) + |0,iV„,ax) I 5 )). (54) 

The recursion algorithm restoring \tjjf) to the vacuum state 
|0, 0 ) 15 ) is schematically shown in Fig. [3] for the maxi¬ 
mum photon number A^max = 2. In Fig. [3ja), we can 
find that all the populations in the Hilbert space spanned by 
{| 0 , 2 )|p), | 2 , 0 )|p)} can be transferred to the state | 1 , 0 )|e) by 
consecutively using transitions “ 01 ”, “ 10 ”, “ 01 ”, and “ 10 ”, 
i.e.. 


|0,2) |5)^|0,l)|e)^|l,l)|g)^ 

|l,0)|e)^ | 2 , 0 )| 3 ). (55) 

After this procedure, as schematically shown in Fig. Ob), all 
the populations on the state 11 , 0 ) | e) can be transferred to the 


state | 0 , 0 )|( 7 ) by consecutively using transitions “00”, “10”, 
and “00”, i.e., 

|l,0)|e)^|l,0)|p)^ 

/-4 f-5 

|0,0)|e)^|0,0)|p). (56) 

/-6 

The total step number is thus / = 7 for generating the NOON 
state (|0,2) -f |2,0 ))/V2 . 

More generally, given an arbitrary A^max, the total step 
number for generating the NOON state in Eq. (l54l) is 

/ = 4At„,ax-l. (57) 

The step number for generating NOON state has been greatly 
reduced in comparison with that for generating an arbitrary 
state in Eq. (IA16b . Obviously, the NOON state can be gen¬ 
erated without using the “11” transition. If we assume the 
Lamb-Dicke parameter gi < 1, which is usually the case even 
in the ultrastrong regime in superconducting circuit QED sys¬ 
tems ll36l - [^ . From Eq. (l22l i. we have the Rabi frequencies 

« P 1 P 2 , OC r]i, and 

J OC 772 - Thus, the transition “11” generally takes 
more time among the four types of transitions employed by 
us. Therefore, our algorithm may show a better efficiency 
for generating NOON sates than generating arbitrary entan¬ 
gled states. This is especially true when the maximal photon 
number A^max is higher and the Lamb-Dicke parameter 77 / is 
smaller. 




FIG. 3: (Color online) Application of the general algorithm to 
generating the NOON state. The notations are the same as those 
in Fig. 1^ (a) Schematic diagram for transferring the population 

in the space {|0, 2)|(/), |2, 0)|(;)} to the state |l,0)|e). This is 
achieved by consecutively using “01”, “10”, “01”, and “10” tran¬ 
sitions. (b) Schematic diagram for transferring the population on the 
state |1, 0)|e) to the state |0,0)|(?). This is achieved by consecutively 
using “00”, “10”, and “00” transitions. 


IV. MINIMIZING THE EFFECT OF UNWANTED TERMS 
A. Theoretical analysis 

In all of the above studies, we make an approximation that 
all unwanted terms have been neglected. However, these ne- 
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glected terms will affect the fidelity of the prepared target 
state. Let us now discuss how to minimize the effect of these 
unwanted terms in Eq. (l20l) on the target state by choosing 
appropriate parameters. In principle, the effects of these un¬ 
wanted terms can be perfectly removed by pulse calibration 
techniques. Here, we study a method to minimize the effect 
of these unwanted terms by choosing the parameters when the 
pulse calibration cannot be used. 

In our algorithm, we have used four interactions 

and i/™, all of them are constructed by the terms with 
the Bessel function Jj (x) in Eq. (l20l) . Here, in the sub¬ 
script of the Bessel function, we also use N to denote —N 
if N > 0. We hope to suppress all the terms with the 
Bessel functions J^i {x) for N' ^ —1. We focus on the 
case a; = 2H/a) ~ 1 considering possible experimental con¬ 
ditions. In this case, only lower order Bessel functions Jq (x), 
J±i (x), and J ±2 (a;) play significant roles. Thus, we need 
only to find proper parameters such that the effect of the terms 
with Jq (x), Ji (x), and J ±2 (x) are negligibly small. Our 
idea is to make those terms nonresonant by properly choos¬ 
ing the parameters and uJz of the qubit, and frequencies wi 
and W 2 of two microwave modes. That is, we assume that the 
frequency of the Zth cavity mode satisfies 

— h (58) 


The conditions in Eqs. (I62l i- (l65l l can be summarized as 


^ k UJgcd- 


( 66 ) 


We can also assume that the longitudinal frequency of the 
qubit is 


U}z = {p + r) Wgcd, 


(67) 


where p is the integer part, and r is the fraction part. To meet 
Eq. dbbl l. there should be 


(123451 

\ ’ 6’ 6’ 6’ 6’ 6 j ■ 


( 68 ) 


The nonresonant terms with Ji, Jq, and J ±2 still have effect 
on the desired time evolution. These effects can be further 
eliminated by decreasing the stark shifts caused by the terms 
with in Eq. (l22l i. The ideal case is 


where k is a positive integer and Wgcd is the greatest common 
divisor of oji and a; 2 . Assuming that the “fcifc 2 ”transition is 
switched on, i.e., the transition detuning = 0, then from 

Eq. ( |23] ). the frequency w of the driving field must satisfy the 
condition 






2 


< 


lnin2 ’ 


(69) 


UJ = UIz + kiLOi + k2U!2- 


(59) or equivalently. 


Erom Eqs. (|2^ . (l58l l. and ( |59] |. the detuning of the term with 
N', k[, k '2 is then given by 

2 

= (TV' + 1 ) 0 ;, + ^ {N'h + k[) hujg,A. (60) 


^ N'n'-^n2 




k\k 2 

ln\n 2 


< |A 


h' h' , 

N' I 


(70) 


Thus the terms with the Bessel function Ji {x) will have the 
detuning 


2 

A]^^ ^ = 2wz -b 'y ^ (ki -\- k'l) huigcd- (61) 

1=1 

We expect that the terms with Ji (x) are nonresonant. Thus, 

k' k' 

the relation that Aj^^ ^ ^ 0 must hold. A simple but sufficient 
condition is 

2^2 kujgcd} (62) 

where k is an integer. Similarly, for the terms with Jq (x), 
J 2 (x), and J 2 (x), the sufficient conditions can be given by 

(63) 

(64) 

(65) 


for N' = 0, 1, ±2, n[ + < Amax, and ni -b 712 < Amax, 

where the constraint condition for ni and nj denotes the work¬ 
ing space of our algorithm. Equation ( |69] | means that the stark 
shifts should be negligibly smaller than the Rabi frequencies 
for state generation. Considering that Amax is the maximum 
photon number of the target state, and using Eq. dbOl l and 
Eq. (|67l), we can obtain 


A 

A 

A 

A 


fcl ^2 
0 

fci ^2 
1 

fel ^2 
2 

fci k2 
2 


> rojgcd or (1 - r) Wgcd 

> (2r - [2rJ)a;gcd or (|■2r■] - 2r-) Wgcd, 

> (3r - [3rJ)a;gcd or (|■3r■] - 3r) Wgcd, 

> rtdgcd or (1 - r) Wgcd- 


(71) 

(72) 

(73) 

(74) 


Here, [xj means x rounded down and [x] means x rounded 


lOz k CCgcd; 

3 dJz k CCgcd; 

Uz k edged- 
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up. We thus reduce Eq. (l70l) to 


B. Numerical simulations 


2 

- < min{r, 1 - rjwgcd, (75) 

2 

- <c niin{2r - [2rJ , |'2r] - 2r}a;gcd, (76) 
2 

- < min{3r - [3rJ , [Sr] - Srjwgcd, (77) 
2 

- < min{r, 1 - rjwgcd, (78) 


a condition much stronger than Eq. ( fTOl i. If Eqs. dTSll-dTSll are 
fulfilled, the nonresonant terms can in principle be suppressed. 
We know from Eq. (|2^ that Eqs. (I75l l- (l78l l can be satisfied 
if, for example, the parameter uj^ of the qubit is tuned suffi¬ 
ciently small, assuming that the reduced driving frequency x 
and Lamb-Dicke parameters rji have been appropriately cho¬ 
sen. 

Beside the terms with Bessel functions Jiy> (x) where N' ^ 
— 1, there are also unwanted terms with the Bessel function 
Jj (x), which, however, also satisfy the resonant condition 

2 

=^(fci'-fci)Z;tdgcd = 0. (79) 

i=\ 

Here, we have used Eqs. ( |59] | and (l60l l to obtain Eq. (|79] |. The 
Lamb-Dicke parameters satisfy the condition r]i = 2gi /l oi 5, 
1 for circuit QED systems even in the ultrastrong regime 113^ 
1^ . Erom Eq. (fSST l. we know that li and I 2 are coprime num¬ 
bers. We can further make li (or I 2 ) sufficiently large. Thus 
the unwanted resonant terms will possess large \k[\ (or 1 ^ 2 !)■ 
In this way, the effects of these terms will be suppressed due 
to the exponential decrease via the term in Eq. (l22l i. The 
condition, that the term Jj (x) is negligibly small, can be sum¬ 
marized as that 1 1 and I 2 should satisfy 

Viri2" < 1. (80) 

We now summarize the condition that minimizes the effects 
of unwanted terms. The parameter of the qubit should sat¬ 
isfy Eq. (l67l i and Eq. (l68l l. However, the parameter of the 
qubit is mainly constrained by current experiments. Eor ex¬ 
ample, typical values of are in the range 1 ~ 5 GHz. 

The frequencies of the cavity modes uji should satisfy Eq. (ISFl l 
and Eq. (I 8 OI 1 . The values of the reduced driving frequency 
X = 2H/u] and Lamb-Dick parameter Pi = 2gi/oji should 
satisfy Eq. ( fTOl i or stronger conditions Eqs. (f75]l-(l78]l. Ap¬ 
propriate values of x and p; can be obtained via numerical 
simulations, which will be discussed below in Sec. lIV Bl 


^Qn'-^n 2 


lnin 2 

ln^n 2 


lnin 2 

Qfcjfcj 

^ 2 n'^n 2 


^fclfe2 

lnin 2 

2 n'^n 2 


^k^k 2 

IniriQ 


We now further numerically simulate the effect of the un¬ 
wanted terms on the generation of target states by using ex¬ 
amples of generating the following two target states 

V'e) = ^ I5) , ( 81 ) 

m+n 2<2 

^n) = ^(|0,2) + |2,0))|5), (82) 

for some given parameters. It is obvious that \ip^) is an en¬ 
tangled state where every state component is evenly occupied. 
We thus call |^e) the evenly-populated state. The state |^n) 
is a two-photon NOON state ll. Both |'0 e) and |^n) pos¬ 
sess a maximum photon number iVmax = 2. The fidelities for 
generating these two states \iIje) and are defined as 


-T^e = 

(^E 


J'n = 

(^N 



(83) 

(84) 


Here, and are respectively the actually generated 
states via the total Hamiltonian in Eq. (|6]). 

We now determine the detailed experimental parameters. 
Erom Eq. dbTl i and Eq. (l68l l. we set r = 3/4, p = 9, and 
Wgcd/27r = 2 GHz, which corresponds to uJz/2 = 19.5 GHz. 
Erom Eq. ( iTQl i or Eqs. (I75ll-(l78ll. the parameter uJxI2.'k should 
be made smaller, e.g., we set uJxI2.'k = 1.2 GHz. The fre¬ 
quency of the kh cavity, i.e., w;, is determined by Eq. (ISST i 
and (IHOl l. Since the microwave fields are usually of several 
gigahertz, here we set 1 1 = 3, and I 2 = 4, thus yielding 
a;i/27r = /iu;gcd/27r = 6 GHz and u;2/27r = /2Wgcd/27r = 8 
GHz. The Lamb-Dicke parameters for the first and second 
cavity modes are set to be identical, i.e.. 


111 = 112 = V- (85) 

We vary the Lamb-Dick parameter p and the reduced driv¬ 
ing frequency x = 2 il/dj to simulate the effect of the un¬ 
wanted terms on the fidelity of the expected target states in 
Eqs. dSTT i and (|82] |. The pulses are taken according to the cal¬ 
culation of Chapter III B. That is, in the i/th step, we use a 
sinusoidal driving with the driving frequency oj^. Since the 
sinusoidal driving lasts for a duration fi,, the driving field can 
be considered as square-windowed sinusoidal signal and thus, 
strictly speaking, is not delta-shaped in the spectrum. The 
simulation results for generating target states in Eq. dSTl i and 
Eq. (l82l i are listed in TableUand Tablellll respectively. We can 
easily find that larger reduced driving strengths x and Lamb- 
Dick parameters p can usually make the fidelity higher. For 
the evenly-populated state |V'e), the largest fidelity 0.939 can 
be obtained at x = 1.7571 and p = 0.3714. However, for the 
NOON state the largest fidelity 0.92 can be obtained at 
X = 2 and p = 0.5429. 
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TABLE I: The fidelities J^e = K'*^b I'*/’e}| of the target state It^e) = 
Eni+n2<2(l/^)l^2,n2)|g) are listed for different values of the 
reduced driving frequency x = and the Lamb-Dicke param¬ 

eter rj — 2gi/ijJi = 2g2lijJ2- Here |'0 b) is the actually gener¬ 
ated state using the total Hamiltonian. We have chosen the longi¬ 
tudinal frequency of the qubit uj^/27r = 19.5 GHz, the transverse 
frequency of the qubit uixl2n = 1.2 GHz, the frequency of the 
first mode cji /27r = 6 GHz and the frequency of the second mode 
aj2/27r = 8 GHz. 



V 

0.2 

0.3714 0.4571 

0.5429 

0.6286 0.7143 


0.3 

0.115 

0.393 

0.519 

0.608 

0.739 

0.675 


0.7857 

0.589 

0.817 

0.872 

0.851 

0.911 

0.852 

X 

1.0286 

0.615 

0.85 

0.876 

0.906 

0.886 

0.857 


1.2714 

0.724 

0.872 

0.886 

0.906 

0.852 

0.878 


1.7571 

0.821 

0.939 

0.899 

0.859 

0.825 

0.837 


2 

0.867 

0.915 

0.838 

0.876 

0.887 

0.859 


TABLE II: The fidelities = |(^n|^n)| of the target state 
I^^noon) = (l/y2)(|0,2}|g} + |2,0)|g)) are listed for different 
values of the reduced driving frequency x = 2H/iij and the Lamb- 
Dicke parameter rj = 2g\/uii = 2y2/uJ2. Here |'0 g} is the actu¬ 
ally generated state using the total Hamiltonian. We have chosen the 
same parameters as in Table. |I] 



V 

0.2 

0.3714 0.4571 

0.5429 

0.6286 0.7143 


0.3 

0.108 

0.34 

0.403 

0.395 

0.461 

0.687 


0.7857 

0.675 

0.815 

0.833 

0.862 

0.829 

0.868 

X 

1.0286 

0.78 

0.857 

0.846 

0.867 

0.883 

0.813 


1.5143 

0.871 

0.877 

0.873 

0.877 

0.787 

0.806 


1.7571 

0.844 

0.918 

0.889 

0.902 

0.862 

0.819 


2 

0.876 

0.909 

0.832 

0.92 

0.853 

0.806 


V. ENVIRONMENTAL EFFECT ON TARGET STATES 

In the above, we only discuss the effect of unwanted terms 
on the generation of target states. We now study the effect of 
dissipation on the fidelities of target states by numerical sim¬ 
ulation for given parameters. When the environmental effect 
is included, the dynamical evolution of the SQC can be de¬ 
scribed by the master equation 


qubit are given by the eigenstates of Eq. (|2]i. If we define 

= (87) 

with V = g,e and g = g,e, and also define 

= W) (mI : (88) 

with iy = g,e and fl = g,e. We can easily verify 

= Ry {0) (Tuy.Rl {0) , 

where Ry {9) = exp [—i6(jyl2), and 9 = arctan(a;a;/wz). In 
Eq. (l86l l. 7 eg is the pure-relaxation rate from the qubit excited 
state to the ground state. Besides, 7 gg and 7 ee are the pure- 
dephasing rates originating from disturbed qubit eigenstates. 
The decay rates of the first and the second cavity fields are 
denoted by ki and K 2 , respectively. 

Using parameters in Sec. lIV Bl and taking the reduced driv¬ 
ing strength x = 2U/w = 1.7571 and Lamb-Dicke param¬ 
eter rj = 0.3714 from Tableland Table HU we find that the 
highest fidelity = 0.939 is achieved for generating the 
evenly-populated state |V’e) in Eq. (ISTT i. and a high fidelity 
= 0.918 is also reached for generating the NOON state 
|-!^n) in Eq. (|82] |. 

We now assume that the decay rates in Eq. (l8^ are taken 
as 7gg/27r = 0, 7ee/27r = 2 MHz, and 7eg/27r = Ki/27r = 
K2/27r = 1 MHz. We assume that the density operators pj| 
and are the actually generated states for the target states 
|'0 e) and IV'n)- Then the fidelities can be redefined as 




pi 


Ip' 


E , 


Rn = \/(V'n 


pi 


V'n)- 


(89) 

(90) 


We perform numerical simulations using the above parame¬ 
ters and obtain = 0.911 and = 0.863. The total time 
for generating {ip's) is Te = 8.9561 ns and that for generating 
I^^n) is 7 n = 10.4451 ns. Both Te and Tn are too small to 
induce significant decoherence at the decay rates specified by 
us. Thus, the fidelity losses induced by dissipation are fairly 
small, which are J^e — — 0 028 for the evenly-populated 

state |'0 e) and — -^n = 0.055 for the NOON state \ip-N)- 


VI. DISCUSSIONS 


P = -i[H,p]+'D [V^O'se] P + T> [v/ET^ee] P 

+ ^ [Vl^^gg] P + R [v/m“i] P + R [y/^a2] p, (86) 

where p and H are the reduced density operator and the 
Hamiltonian of the whole system, respectively. The total 
Hamiltonian has been given in Eq. (| 6 l). The compact notation 
T>[c] p = {2cpc^ — c^cp — pc^c) /2 represents the Lindblad- 
type dissipation. We have noted that {I 5 ), |e)} is the basis of 
az, but the qubit dissipation is determined by the qubit ba¬ 
sis Up), |e)}. The ground (|p)) and excited (|e)) states of the 


We now discuss the advantages and disadvantages between 
our methods and the previous ones ItoI - It^ for generating ar¬ 
bitrarily entangled states of two microwave fields or two vi¬ 
brational modes. 

The brief comparison between these methods is listed in Ta¬ 
ble. ini] In detail. Ref. provided an algorithm to generate 
arbitrarily entangled states of two vibrational modes. But due 
to population leakage outside the original space, it takes an 
exponential complexity of the number of steps. The succeed¬ 
ing proposals llEll - l^ overcome the exponential drawback in 
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TABLE III: Comparison of different methods for generating arbitrar¬ 
ily entangled states of two-mode bosonic fields. We use Pop. Leak., 
St. Eff. and Mult. Proc. to denote population leakage, the Stark 
effect, multiboson processes, respectively. No. At. Lev. is used to 
denote the number of atomic energy levels. For example, 2 denotes 
two energy levels when the state is generated. We use ”Yes” or ”No” 
to denote whether the population leakage (Stark effect) occurs (are 
used) or not. Meanwhile, L (or H). Pn. No. means multiphonon pro¬ 
cesses of low (or high) phonon number, however, L (or H). Pt. No. 
means multiphoton processes of low (or high) photon number. 



Pop. Leak. 

No. At. Lev. 

St. Eff. 

Mult. Proc. 

Ref. [20] 

Yes 

2 

No 

L. Pn. No. 

Ref. \2V\ 

No 

3 

No 

L. Pn. No. 

Ref. [2^ 

No 

3 

No 

H. Pn. No. 

Ref. [23J 

No 

2 

Yes 

None 

Ref. [2T] 

No 

2 

No 

H. Pn. No. 

Refs. [25-27] 

No 

2 

Yes 

None 

Our proposal 

No 

2 

No 

L. Pt. No. 


several ways: (1) A third atomic level is used to shield oscilla¬ 
tions that cause population leakage ll^l^ . But the disadvan¬ 
tage is that higher energy levels of systems usually have larger 
decay rates, which inevitably reduce the fidelities of the tar¬ 
get states. (2) Boson-number-dependent Stark effects are used 
to realize independent operations of particular states l25l - 
IH. But the disadvantage is that the detunings of nonresonant 
terms are usually less by one order of the coupling strengths 
between the two-level system and boson modes. This means 
that the Rabi frequencies are smaller, and the longer genera¬ 
tion time is required. (3) Multiphoton processes of high pho¬ 
ton number are used to shield oscillations that cause popu¬ 
lation leakage or reduce the number of steps I^ . But 
the disadvantage is that if the coupling strengths between the 
atom and cavity fields are not high enough, then the Rabi fre¬ 
quencies become small, especially for states with high photon 
numbers, which obviously indicates longer generation time. 

Besides the advantage that there is no population leakage, 
our method has also the following advantages compared with 
previous ones. ll2l1 - l^ : (1) It only uses the two energy lev¬ 
els of the qubit. Thus, the fidelities of the target states should 
be higher because there is no other auxiliary energy levels. 
(2) The detunings of the nonresonant terms are in the order 
of the resonator frequencies. They are usually bigger than 
the coupling strengths between the qubit and resonator modes. 
Thus the Rabi frequency can be made bigger than those using 
boson-number-dependent Stark effects. (3) We use multipho¬ 
ton processes of low photon number, i.e., one photon at most 
is converted for either mode. Thus the Rabi frequency can be 
bigger than those using multiphoton processes of higher pho¬ 
ton number,especially when the coupling strengths between 
the qubit and cavity modes are not very big. Of course, 
stronger couplings will further enhance the Rabi frequencies 
and hence reduce the generation time. 

We point out that the real supercoducting qubit circuits are 
mutilevel systems, the information leakage to higher levels is 
not avoidable. However, the leakage can be neglected when 


the transition frequency between the first excited state and the 
second excited state is much larger than the qubit frequency. 
For example, in the flux qubit circuits, due to its large anhar- 
monicity of energy levels, the information leakage is negli¬ 
gibly small. However, for the transmon and phase qubit, the 
anharmonicity is very weak. Thus, the pulse should be care¬ 
fully calibrated to avoid information leakage to higher levels. 
The pulse calibration can be done as in Ref. 1^ . 

We now compare the differences between our algorithm 
and other ones for generating NOON states. Ref. IMBI] uses 
mutliphoton processes to generate NOON states. In supercon¬ 
ducting systems, this means a low generation efficiency if the 
Lamb-Dicke parameter is not sufficiently big. Ref. 11^ uses 
synchronization technology to generate NOON states, but the 
time duration for synchronization between two steps can be 
quite long and there exists inevitably information leakage. 
Ref. and its experimental realization ll43l] use two phase 
qubits with three active energy levels to generate NOON states 
of two cavity modes. The experimental setup is complex and 
the high energy levels of qubits will reduce the decoherence 
time. Ref. uses photon number-dependent Stark effects 
to achieve independent operations. Thus the Rabi frequency 
is smaller than the qubit-cavity coupling strengths. Ref. 114^ 
requires that two qubits be initially prepared in a Bell state 
and Anally get decoupled from the qubits and cavity fields. 
Ref. ii^ll uses one qubit but still needs one additional level 
to shield unwanted resonances. More recently. Ref. uses 
one qubit of four levels which resonantly interacts with two 
resonators simultaneously to speed up the generation process 
of NOON states. 

When applied to generating NOON states, our algorithm 
has new features besides the common advantages for gener¬ 
ating arbitrary two-mode photon states: (1) Only carrier pro¬ 
cesses and one-photon processes are used. In this case, 
even though the coupling strengths between the qubit and cav¬ 
ity modes are small, large Rabi frequencies can still be ob¬ 
tained. (2) The number of steps is reduced to linear depen¬ 
dence on the maximum photon number. These advantages 
indicate less generation time and thus guarantee a higher effi¬ 
ciency than preceding methods. 

Now we discuss the experimental feasibility of our scheme. 
TableUand Tablellllshow that without pulse calibration, higher 
fidelities can be achieved at bigger Lamb-Dicke parameters 77 
and reduced driving frequencies x. These values are already in 
the ultrastrong regime. Ref. has reported ultrastrong cou¬ 
plings between three resonator modes and a flux qubit, where 
the Lamb-Dicke parameter 77 can reach as high as 0.236. In the 
ultrastrong regime, Rabi frequencies can be made to approach 
the magnitude of uj^, which usually ranges from 1 to 5 GHz. 
The decay rates of the qubit and cavity fields are usually in the 
magnitude of megahertz. Thus the dissipation has small effect 
on the fidelities of target states. For singe-mode microwave 
fields, Fock states with up to six photons lH] and Fock state 
superpositions d have been experimentally demonstrated 
using phase qubits. The NOON state up to 3 photons has also 
been experimentally reported [|43|. We thus hope that our 
proposal is also experimentally feasible in the near future. 
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VII. CONCLUSIONS 


denoted as 


In summary, we have proposed an approach to generate 
arbitrary superpositions of photon states of two microwave 
fields in two separated cavities. Our method mainly depends 
on the coexistence of transverse and longitudinal couplings 
between the qubit and cavity fields. Employing the longitudi¬ 
nal couplings, we derive a Hamiltonian which is similar to that 
of trapped ions interacting with two vibrational modes 0. 
Using four simple interaction Hamiltonians derived from the 
longitudinal coupling, we design the state generation algo¬ 
rithm. Ou r alg orithm can be regarded as the improved version 
of that in when the transverse and longitudinal couplings 
coexist in circuit QED systems. But it has remedied the draw¬ 
back that the number of steps exponentially depends on the 
maximal photon number, which is replaced by a quadratic de¬ 
pendence. Compared with previous ones with quadratic com¬ 
plexity, our algorithm does not require atomic energy levels 
higher than two ll^ |2^ . boson-number-dependent stark ef¬ 
fects or multiboson processes of high boson num- 

When applied to the generation of NOON states, whose 
engineering has been extensively studied lE^-Eitl, our algo¬ 
rithm needs only carrier and one-photon sideband transitions. 
Meanwhile, the number of steps only linearly depends on the 
maximum photon numbers. In fact, these properties for gen¬ 
erating NOON states can be generalized to any states with a 
constant total photon number of both modes. 

We have also discussed how to avoid the effect of unwanted 
terms on the generation of target state. Our numerical results 
show that fidelities above 0.91 can be reached in the ultra¬ 
strong regime for the two-photon evenly-populated state and 
NOON state when the environmental effect is neglected. The 
generation time can be very short, in which case, the envi¬ 
ronment has small effect on fidelities of the target states. We 
here note that due to the similarity of two-mode interaction 
Hamiltonians, the alg orithm using two-mode multi-phonon 
processes in Ref. can be directly applied into our model. 
Thus, two-mode Eock states with high photon numbers can be 
generated with just two steps as one-mode Eock states in 0. 

We have noted that our method for generating NOON 
states is similar to a recent algorithm simplified from the one 
which employs Stark effects to generate arbitrary entangled 
states ll2^ . 
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Appendix A: Universal algorithm for arbitrary maximum 
photon numbers 


'^ 2 Arn,ax = {\ni,n 2 ) \g) \ni + n 2 < iV^ax} • (Al) 

We will implement 27Vniax procedures, with each procedure 
containing some steps. 

In the 1st procedure, we aim to clear the populations in the 
subspace 

= { 1 ^ 11 ^ 2 ) |g) +n2 = A^max} (A2) 

of ’H2Arn,ax ■ This can be achieved via alternatively switching 
A^max “10” transitions and Wmax “01” transitions, i.e.. 


|0,lV„.ax)|5) ^ |0,lV„.ax-l)|e) |l,7V„,ax-l)|<?)' 

> lA^max - 1, 0) |e) - |Wn,ax, 0) \g) . 


/-2Afmax-|-2 


/-2Afmax-|-l 


(A3) 


to transfer the populations in the subspace ^2Arn,ax 

to the state |Wmax — 1,0) |e). Thus, ’H2Arn,ax 1^ reduced to the 
state space 'H2Ar„ax-i where 

'^2Arn,ax-i ={\ni,n2) \g) \ni + n2 < Wmax - 1 } 

U {|ni,n 2 ) |e) \ni +n 2 < A^max - 2} 
U{|Wn,ax-l,0)|e)}. (A4) 


In the 2/ith procedure, we aim to clear the populations in the 
subspace {|Wmax - /i,0) |e)} of 'H 2 Ar„ax- 2 /i-ti-_This can be 
achieved via alternatively switching Wmax—M “11” transitions 
and A^max — g “00” transitions, i.e.. 


|0,Wn,ax-M)l5) |l,Wn,ax-M-l)|e) 

f-N 2 ti 


00 




> |l,A^max-M - 1) Iff)- 


11 


f-N 2 ^- 2 {Nm^^-fl )-\-2 
00 


/ — A^2— 2 ( A^m ax — M ) +1 


> Amax - /i,0) |e) 

> Amax - /i,0) Ip) . 


(A5) 


with A^2 m = 2Wniax + (4A^max - 2/r - 1) (p - 1), to transfer 
the populations in the subspace 

^2Ar„„-2M-ti ={\ni,n2) \g) |ni -f n2 = g.,n2 ^ 0 } 

U{|W„,ax-/i,0)|e)} (A6) 

of 'H 2 Nr^^^ - 2 n+i to the state IW^ax — ft,0) |p). Here Thus, 
'^2Af„ax-2M-l-i is reduced to the state space 'H2Af„,ax-2At where 

’^2Afn,„-2M ={\ni,n2) Ip) \ni + n2 < A^^ax - M - 1} 

U {|ni,n 2 ) |e) \ni + n 2 < A^^ax - p - 1} 
U{|lV,„ax-p,0)|p)}. (A7) 


The basic principle of our algorithm is state space reduc¬ 
tion. The state space of the target state in Eq. (l33]) can be 


In the (2/i + l)th procedure, we aim to clear the popu¬ 
lations in the subspace {lA^^ax - P, 0) |p)} of 'H 2 Af„ax- 2 /i- 
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This can be achieved via alternatively switching -/Vmax — M of "Hi. This can be achieved via switching one “00” transi- 
“10” transitions and A^max — /r — 1 “01” transitions, i.e., tions, i.e., 


|0,iVmax - M - 1) |e) 


10 


/-W2, 


|li A^max 1) |<?) 


|0,0)|e) ^ |0,0)|g). (A12) 


01 




> |l,iVniax - - 2) |e) • 


10 


/ —1 — 2 ( A/'max — m)+3 
01 


/ —A^2/2 + 1 — 2(A^inax—M) + 2 


> |A^max - At - 1,0) |e) 
|Atniax ~ Atj 0) li?) ■ 


(A8) 


withA’2^+1 = At2^+2(Atniax —m)’ to transfer the populations 
in the subspace 

'^2Ar„ax-2M ={\ni,n2) \e) \ni +712 = At^ax - ti-l,n2^0} 
U{|iVa,ax-/i,0)|p)} (A9) 

of T-L2Nm^^ -2fj, to the state |At„iax — — ^, 0 ) \e). Thus, 

'^2Af„,ax-2M is reduced to the state space 'H2Af„ax-2M-i where 

'^2Arn,ax-2M-l ={|«1,«2) \g) \ni + 712 < A^max “ A^ “ 1 } 

U {|ni,n 2 ) |e) |ni + 112 < At^ax - /i - 2} 

U{|At„,ax-At-l,0 )|e)}. (AlO) 


to transfer the populations in the subspace H'l of TLi to the 
state |0,0) Ip). Thus, TLi is reduced to the state space T-Lq 
where 


^0 = {|0,0)|p)}, (A13) 

is namely the intial state space for j-f/jo) in Eq. (l34l) . 

Base on the above discussion, we now calculate the number 
of steps to generate the target state in Eq. (l3^ with an arbitrary 
maximum photon number A^max- For alternatively switching 
“11” and “00” transitions, we need f^i qq steps, given by 

^max 1 

/iLoo = l+ E 2N = Nl^^-N^,^ + l. (A14) 

iV=l 


For alternatively switching “10” and “01” transitions, we need 
/io,oi steps, given by 

^max 

/io,oI = l+ E (2At-l)=lVLx + l- (A15) 

Ar=l 


In the (2 A^max)th procedure, we aim to clear the populations 
in the subspace 


Therefore, the total step number / is 

/ = /ll,00 + /l0,0l = ~ A^max + 2. 


H;={|0,0)|e)} 


(A11) for generating the target state in Eq. (l3^ . 


(Alb) 


[1] A. Blais, R. S. Huang, A. Wallraff, S. M. Girvin, and R. J. 
Schoelkopf, Phys. Rev. A 69 , 062320 (2004). 

[2] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R. S. Huang, 
J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 
(London) 431 , 162 (2004). 

[3] S. Ritter, C. Ndlleke, C. Hahn, A. Reiserer, A. Neuzner, M. 
Uphoff, M. Miicke, E. Figueroa, J. Bochmann, and G. Rempe, 
Nature (London) 484 , 195 (2012). 

[4] H. J. Kimble, Nature (London) 453, 1023 (2008). 

[5] M. A. Nielsen and 1. L. Chuang, Quantum Computation and 
Quantum Information (Cambridge university press, 2010). 

[6] M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge 
university press, 1997). 

[7] Y. X. Liu, L. F. Wei, and F. Nori, Europhys. Lett. 67, 941 
(2014). 

[8] Yu-xi Liu, L. F. Wei, and F. Nori, Phys. Rev. A 71, 063820 
(2005). 

[9] M. Hofheinz, E. M. Weig, M. Ansmann, R. C. Bialczak, E. 
Lucero, M. Neeley, H.Wang, J.M.Martinis, and A. N. Cleland, 
Nature (London) 454 , 310 (2008). 

[10] M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, E. 
Lucero, M. Neeley, A. D. O’Connell, D. Sank, J. Wenner, J. M. 
Martinis, and A. N. Cleland, Nature (London) 459 , 546 (2009). 

[11] K. Vogel, V. M. Akulin, and W. P. Schleich, Phys. Rev. Lett. 71, 


1816 (1993). 

[12] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73 , 
357 (2001). 

[13] J. Q. You and F. Nori, Phys. Today 58, 42 (2005). 

[14] G. Wendin and V. S. Shumeiko, Low Temp. Phys. 33, 724 
(2007). 

[15] R. Schoelkopf and S. Girvin, Nature (London) 451 , 664 (2008). 

[16] J. Clarke and F. K. Wilhelm, Nature (London) 453, 1031 
(2008). 

[17] J. Q. You and F. Nori, Nature (London) 474 , 589 (2011). 

[18] 1. Buluta, S. Ashhab, and F. Nori, Rep. Prog. Phys. 74, 104401 

( 2011 ). 

[19] Z. L. Xiang, S. Ashhab, J. Q. You, and F. Nori, Rev. Mod. Phys. 
85, 623 (2013). 

[20] S. A. Gardiner, J. 1. Cirac, and P. Zoller, Phys. Rev. A 55, 1683 
(1997). 

[21] G. Drobny, B. Hladky, and V. Buzek, Phys. Rev. A 58, 2481 
(1998). 

[22] S. B. Zheng, Phys. Rev. A 63, 015801 (2000). 

[23] B. Kneer and C. K. Law, Phys. Rev. A 57, 2096 (1998). 

[24] X. B. Zou, K. Pahlke, and W. Mathis, Phys. Rev. A 65 , 045801 
( 2002 ). 

[25] F. W. Strauch, K. Jacobs, and R. W. Simmonds, Phys. Rev. 
Lett. 105 , 050501 (2010). 







14 


[26] F. W. Strauch, D. Onyango, K. Jacobs, and R. W. Simmonds, 
Phys. Rev. A 85, 022335 (2012). 

[27] R. Sharma and F. W. Strauch, arXiv: 1503.02157 (2015). 

[28] Y. X. Liu, J. Q. You, L. F. Wei, C. P. Sun, and F. Nori, Phys. 
Rev. Lett. 95, 087001 (2005). 

[29] Y. X. Liu, C. X. Yang, H. C. Sun, and X. B. Wang, New J. 
Phys. 16, 015031 (2014). 

[30] Y. J. Zhao, Y. L. Liu, Y. X. Liu, and F. Nori, Phys. Rev. A 91, 
053820 (2015). 

[31] M. Sasura and V. Buzek, J. Mod. Opt. 49, 1593 (2002). 

[32] L. F. Wei, Y. X. Liu, and F. Nori, Phys. Rev. A 70, 063801 
(2004). 

[33] J. Steinbach, J. Twamley, and P. L. Knight, Phys. Rev. A 56, 
4815 (1997). 

[34] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, Phys. 
Rev. Lett. 89, 117901 (2002). 

[35] J. Lisenfeld, A. Lukashenko, M. Ansmann, J. M. Martinis, and 
A. V. Ustinov, Phys. Rev. Lett. 99, 170504 (2007). 

[36] T. Niemczyk, F. Deppe, H. Huebl, E. Menzel, F. Hocke, M. 
Schwarz, J. Garcia-Ripoll, D. Zueco, T. Hummer, and E. 
Solano, Nat. Phys. 6, 772 (2010). 

[37] P. Eorn-Diaz, J. Lisenfeld, D. Marcos, J. J. Garci'a-Ripoll , E. 
Solano, C. J. P. M. Harmans, and J. E. Mooij, Phys. Rev. Lett. 
105, 237001 (2010). 


[38] R. Stassi, A. Ridolfo, O. Di Stefano, M. J. Hartmann, and S. 
Savasta, Phys. Rev. Lett. 110, 243601 (2013). 

[39] M. Steffen, J. M. Martinis, and I. L. Chuang, Phys. Rev. B 68, 
224518 (2003). 

[40] X. B. Zou, J. Kim, and H. W. Lee, Phys. Rev. A 63, 065801 

( 2001 ). 

[41] X. W. Xu, Y. J. Zhao, and Y. X. Liu, Phys. Rev. A 88, 022325 
(2013). 

[42] S. T. Merkel and E. K. Wilhelm, New J. Phys. 12, 093036 

( 2010 ). 

[43] H. Wang, M. Mariantoni, R. C. Bialczak, M. Lenander, E. 
Lucero, M. Neeley, A. D. O’Connell, D. Sank, M. Weides, J. 
Wenner, T. Yamamoto, Y. Yin, J. Zhao, J. M. Martinis, and A. 
N. Cleland, Phys. Rev. Lett. 106, 060401 (2011). 

[44] E. W. Strauch, Phys. Rev. Lett. 109, 210501 (2012). 

[45] Q. P. Su, C. P. Yang, and S. B. Zheng, Sci. Rep. 4, 3898 (2014). 

[46] S.-J. Xiong, Z. Sun, J.-M. Liu, T. Liu, and C.-P. Yang, Opt. Lett. 
40, 2221 (2015). 

[47] C. P. Yang, S. 1. Chu, and S. Han, Phys. Rev. A 67, 042311 
(2003). 

[48] Z. H. Peng, Y. X. Liu, Y. Nakamura, and J. S. Tsai, Phys. Rev. 
B 85, 024537 (2012). 


